The greenhouse effect: the impacts of carbon dioxide (CO2), ultraviolet-B (UV-B) radiation and ozone (O3) on vegetation (crops)

Man's influence on the greenhouse effect, the heating of the atmosphere due to increasing concentrations of tropospheric trace gases, is of much international concern. Among the climatic variables, elevated levels of carbon dioxide (CO₂), ultraviolet-B (UV-B) radiation and ozone (O₃) are known to have a direct effect on vegetation. Our current knowledge of these effects is mainly based on studies involving single stress mode. Thus, the joint effects of CO₂, UV-B and O₃ on vegetation are poorly understood. Nevertheless, based on the literature analysis of plant response to individual stress factors, it can be concluded that sorghum, pea, bean, potato, oat, lettuce, cucumber, rice and tomato are among the crop species potentially sensitive to the joint effects of the aforementioned three variables. Similar information for tree species is essentially lacking.At least with some climatic variables such as O₃, present modeling efforts of cause-effect relationships have proven to be controversial. While at a regional geographic scale ambient CO₂ concentrations appear to be relatively homogeneous, ambient concentrations of O₃ exhibit significant temporal and spatial variability. Because of the protective action of O₃ against UV-B, similar but inverse temporal and spatial variability is expected in the surface levels of UV-B. Thus, future experimental designs should consider these exposure dynamics and modeling cuase-effect relationships should be directed to stochastic processes.     

A data-driven framework to model the organism–environment system

Organisms modify their development and function in response to the environment. At the same time, the environment is modified by the activities of the organism. Despite the ubiquity of such dynamical interactions in nature, it remains challenging to develop models that accurately represent them, and that can be fitted using data. These features are desirable when modeling phenomena such as phenotypic plasticity, to generate quantitative predictions of how the system will respond to environmental signals of different magnitude or at different times, for example, during ontogeny. Here, we explain a modeling framework that represents the organism and environment as a single coupled dynamical system in terms of inputs and outputs. Inputs are external signals, and outputs are measurements of the system in time. The framework uses time-series data of inputs and outputs to fit a nonlinear black-box model that allows to predict how the system will respond to novel input signals. The framework has three key properties: it captures the dynamical nature of the organism–environment system, it can be fitted with data, and it can be applied without detailed knowledge of the system. We study phenotypic plasticity using in silico experiments and demonstrate that the framework predicts the response to novel environmental signals. The framework allows us to model plasticity as a dynamical property that changes in time during ontogeny, reflecting the well-known fact that organisms are more or less plastic at different developmental stages.     

Comparison of conformational characteristics in structurally similar protein pairs.

Although it is known that three-dimensional structure is well conserved during the evolutionary development of proteins, there have been few studies that consider other parameters apart from divergence of the main-chain coordinates. In this study, we align the structures of 90 pairs of homologous proteins having sequence identities ranging from 5 to 100%. Their structures are compared as a function of sequence identity, including not only consideration of C alpha coordinates but also accessibility, Ooi numbers, secondary structure, and side-chain angles. We discuss how these properties change as the sequences become less similar. This will be of practical use in homology modeling, especially for modeling very distantly related or analogous proteins. We also consider how the average size and number of insertions and deletions vary as sequences diverge. This study presents further quantitative evidence that structure is remarkably well conserved in detail, as well as at the topological level, even when the sequences do not show similarity that is significant statistically.     

An automated method for modeling proteins on known templates using distance geometry.

We present an automated method incorporated into a software package, FOLDER, to fold a protein sequence on a given three-dimensional (3D) template. Starting with the sequence alignment of a family of homologous proteins, tertiary structures are modeled using the known 3D structure of one member of the family as a template. Homologous interatomic distances from the template are used as constraints. For nonhomologous regions in the model protein, the lower and the upper bounds for the interatomic distances are imposed by steric constraints and the globular dimensions of the template, respectively. Distance geometry is used to embed an ensemble of structures consistent with these distance bounds. Structures are selected from this ensemble based on minimal distance error criteria, after a penalty function optimization step. These structures are then refined using energy optimization methods. The method is tested by simulating the alpha-chain of horse hemoglobin using the alpha-chain of human hemoglobin as the template and by comparing the generated models with the crystal structure of the alpha-chain of horse hemoglobin. We also test the packing efficiency of this method by reconstructing the atomic positions of the interior side chains beyond C beta atoms of a protein domain from a known 3D structure. In both test cases, models retain the template constraints and any additionally imposed constraints while the packing of the interior residues is optimized with no short contacts or bond deformations. To demonstrate the use of this method in simulating structures of proteins with nonhomologous disulfides, we construct a model of murine interleukin (IL)-4 using the NMR structure of human IL-4 as the template. The resulting geometry of the nonhomologous disulfide in the model structure for murine IL-4 is consistent with standard disulfide geometry.     

Structural analysis based on state-space modeling.

A new method has been developed to compute the probability that each amino acid in a protein sequence is in a particular secondary structural element. Each of these probabilities is computed using the entire sequence and a set of predefined structural class models. This set of structural classes is patterned after Jane Richardson''s taxonomy for the domains of globular proteins. For each structural class considered, a mathematical model is constructed to represent constraints on the pattern of secondary structural elements characteristic of that class. These are stochastic models having discrete state spaces (referred to as hidden Markov models by researchers in signal processing and automatic speech recognition). Each model is a mathematical generator of amino acid sequences; the sequence under consideration is modeled as having been generated by one model in the set of candidates. The probability that each model generated the given sequence is computed using a filtering algorithm. The protein is then classified as belonging to the structural class having the most probable model. The secondary structure of the sequence is then analyzed using a "smoothing" algorithm that is optimal for that structural class model. For each residue position in the sequence, the smoother computes the probability that the residue is contained within each of the defined secondary structural elements of the model. This method has two important advantages: (1) the probability of each residue being in each of the modeled secondary structural elements is computed using the totality of the amino acid sequence, and (2) these probabilities are consistent with prior knowledge of realizable domain folds as encoded in each model. As an example of the method''s utility, we present its application to flavodoxin, a prototypical alpha/beta protein having a central beta-sheet, and to thioredoxin, which belongs to a similar structural class but shares no significant sequence similarity.     

Modeling of protein loops by simulated annealing.

A method is presented to model loops of protein to be used in homology modeling of proteins. This method employs the ESAP program of Higo et al. (Higo, J., Collura, V., & Garnier, J., 1992, Biopolymers 32, 33-43) and is based on a fast Monte Carlo simulation and a simulated annealing algorithm. The method is tested on different loops or peptide segments from immunoglobulin, bovine pancreatic trypsin inhibitor, and bovine trypsin. The predicted structure is obtained from the ensemble average of the coordinates of the Monte Carlo simulation at 300 K, which exhibits the lowest internal energy. The starting conformation of the loop prior to modeling is chosen to be completely extended, and a closing harmonic potential is applied to N, CA, C, and O atoms of the terminal residues. A rigid geometry potential of Robson and Platt (1986, J. Mol. Biol. 188, 259-281) with a united atom representation is used. This we demonstrate to yield a loop structure with good hydrogen bonding and torsion angles in the allowed regions of the Ramachandran map. The average accuracy of the modeling evaluated on the eight modeled loops is 1 A root mean square deviation (rmsd) for the backbone atoms and 2.3 A rmsd for all heavy atoms.     

Surveying nonvisual arrestins reveals allosteric interactions between functional sites

Arrestins are important scaffolding proteins that are expressed in all vertebrate animals. They regulate cell-signaling events upon binding to active G-protein coupled receptors (GPCR) and trigger endocytosis of active GPCRs. While many of the functional sites on arrestins have been characterized, the question of how these sites interact is unanswered. We used anisotropic network modeling (ANM) together with our covariance compliment techniques to survey all the available structures of the nonvisual arrestins to map how structural changes and protein-binding affect their structural dynamics. We found that activation and clathrin binding have a marked effect on arrestin dynamics, and that these dynamics changes are localized to a small number of distant functional sites. These sites include α-helix 1, the lariat loop, nuclear localization domain, and the C-domain β-sheets on the C-loop side. Our techniques suggest that clathrin binding and/or GPCR activation of arrestin perturb the dynamics of these sites independent of structural changes.     

Dynamic genome-scale modeling of Saccharomyces cerevisiae unravels mechanisms for ester formation during alcoholic fermentation

Fermentation employing Saccharomyces cerevisiae has produced alcoholic beverages and bread for millennia. More recently, S. cerevisiae has been used to manufacture specific metabolites for the food, pharmaceutical, and cosmetic industries. Among the most important of these metabolites are compounds associated with desirable aromas and flavors, including higher alcohols and esters. Although the physiology of yeast has been well-studied, its metabolic modulation leading to aroma production in relevant industrial scenarios such as winemaking is still unclear. Here we ask what are the underlying metabolic mechanisms that explain the conserved and varying behavior of different yeasts regarding aroma formation under enological conditions? We employed dynamic flux balance analysis (dFBA) to answer this key question using the latest genome-scale metabolic model (GEM) of S. cerevisiae. The model revealed several conserved mechanisms among wine yeasts, for example, acetate ester formation is dependent on intracellular metabolic acetyl-CoA/CoA levels, and the formation of ethyl esters facilitates the removal of toxic fatty acids from cells using CoA. Species-specific mechanisms were also found, such as a preference for the shikimate pathway leading to more 2-phenylethanol production in the Opale strain as well as strain behavior varying notably during the carbohydrate accumulation phase and carbohydrate accumulation inducing redox restrictions during a later cell growth phase for strain Uvaferm. In conclusion, our new metabolic model of yeast under enological conditions revealed key metabolic mechanisms in wine yeasts, which will aid future research strategies to optimize their behavior in industrial settings.     

Modeling an immune response to influenza A virus infection in alveolar epithelial cells

Influenza A viruses (IAV) have been the cause of several influenza pandemics in history and are a significant threat for the next global pandemic. Hospitalized influenza patients often have excess interferon production and a dysregulated immune response to the IAV infection. Obtaining a better understanding of the mechanisms of IAV infection that induce these harmful effects would help drug developers and health professionals create more effective treatments for IAV infection and improve patient outcomes. IAV stimulates viral sensors and receptors expressed by alveolar epithelial cells, like RIG-I and toll-like receptor 3 (TLR3). These two pathways coordinate with one another to induce expression of type III interferons to combat the infection. Presented here is a queuing theory-based model of these pathways that was designed to analyze the timing and amount of interferons produced in response to IAV single stranded RNA and double-stranded RNA detection. The model accurately represents biological data showing the necessary coordination of the RIG-I and TLR3 pathways for effective interferon production. This model can serve as the framework for future studies of IAV infection and identify new targets for potential treatments.